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EXACT  RESULT  FOR  THE  GRAZING  ANGLE  OF  SPECULAR 
REFLECTION  FROM  A  SPHERE 


1.  Introduction 

When  the  angle  of  incidence  and  the  angle  of  reflection  of  a  wave  or  signal  equal,  we 
say  that  we  have  specular  reflection.  The  angle  of  specular  reflection  or  grazing  angle  from 
a  perfectly  flat  surface  is  easy  to  compute,  but  this  same  angle  from  a  perfect  sphere  is  com¬ 
puted  with  some  difficulty  and  no  exact  formula  for  it  is  presently  known.  The  problem  for 
the  sphere  goes  back  almost  half  a  century  to  W.T.  Fishback  [1]  who  gave  cin  approximation 
for  the  grazing  angle. 

Let  r  be  the  radius  of  a  sphere,  hi  and  h2  respectively  the  perpendicular  heights  of  an 
observer  and  source  above  the  sphere  and  <f>  the  angle  (measured  from  the  center  of  the 
sphere)  between  the  observer  and  source.  If  we  trace  a  wave  from  soiirce  to  observer  (or 
vice  versa),  let  V’  he  the  specular  angle  of  reflection  from  the  sphere  (see  Figure  1).  Miller 
and  Vegh  [2]  gave  an  easy  iterative  method  for  computing  the  grazing  angle  to  any  degree  of 
acairacy.  They  also  showed  that  in  principle  a  formula  for  the  grazing  angle  can  be  foimd 
by  computing  the  roots  of  the  self-inversive  quairtic  polynomial  equation: 

az*  +  I3z^ a  =  0  (1) 


where 


Q  =  e^(e^  -kiki) 

P  =  k^  +  ]4-2k,k2e^ 

7  =  2(fcJ  -I-  -  kik2  cos  <^  -  1)  ; 


and 

hi  =  r/(r -f- hi)  ,  t  =  l,2 
0  <  0  <  COS“^(hl)  +  COS“^(fc2)  <  ’T  • 

If  the  latter  inequality  does  not  hold,  reflection  cannot  occur  and  no  grazing  angle  exists. 
Further,  if  z,  designates  any  of  the  roots  of  equation  (1)  amd 

V'.  =  i  cos“^(Re  2.) 
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they  showed  that  the  uxiique  value  ^  of  V*.  that  satisfies 


4>  +  2^  =  cos  ^(fci  cosV»)  +  co8~^{fca  cos  V') 


(2) 


is  the  grazing  angle. 

Recently,  Secrest  proved  the  conjecture  of  Miller  [3]  that  the  roots  of  equation  (1)  have 
absolute  value  one.  Moreover,  Secrest  showed  that  for 

0  <  fci,/s2  <  1  ,  fci  #  ha  ,  0  <  ^  <  TT  (3) 

four  distinct  roots  of  equation  (1)  enst.  It  was  noted  there  also  that  the  roots  themselves 
are  unknown  in  general. 

In  this  paper,  we  shall  use  the  fact  that  equation  (1)  has  distinct  roots  of  absolute  value 
one  when  the  inequalities  (3)  hold  to  factor  the  quartic  polynomial,  thus  giving  an  explicit 
formula  for  the  grazing  angle.  When  h  =  hi  =  hj  the  grazing  angle  is  given  by 

=  tan"^  (cot  ^/2  —  h  cosec  <f>/2)  . 


2.  Preliminaries 

In  equation  (1)  we  make  the  transformation 

2  =  {a/ay^*8  =  “»“s  (4) 

thereby  rotating  the  roots  thru  an  angle  of  (—1/2)  arga.  Thus  we  arrive  at 

3*  +  +  bs^  +  as  +  1  =  0  (5) 

where 

a  =  /3(|a(a)*’^’ 

*  =  7/\o\  . 

Since  equation  (1)  has  distinct  roots  of  absolute  value  one,  then  so  does  equation  (5)  which 
we  factor  into  two  quadratics: 

(s’  +  +  e’^)(s’  +  U2^  +  e”’^)  =  0  (6) 

where  ui,U2,6  are  to  be  determined.  Writing  equation  (6) 

(e~^s’  +  Uje'^s  +  e^)(c^s’  -f  uac^s  +  e~^)  =  0 
since  each  quadratic  factor  has  zeros  on  the  unit  circle,  we  must  have 

ui  =  me^  ,  U2  =  (7) 


where  /ii,fi2  are  real  and  ||ii|  <  2,  |/i2|  ^  2.  Now  comparing  the  coefficients  of  like  powers 
of  s  in  equations  (5)  and  (6)  and  using  equations  (7)  we  obtain 


a  =  fiie^  +  fj.2e  ^ 

(8) 

b  =  +  H1H2  ■ 

(9) 

Solving  equation  (8)  for  /ii  and  H2  we  deduce 

1  /Re  0  ^  Im  a'\ 

(10) 

2  \cos5  ^  sin 5/ 

4 


(11) 


and  defining 


1  /Re  a  Im  a\ 

2  \cos5  sin^/  ’ 


t  =  cos’ 6  ,  0  <6  <  7r/2 


we  find  from  equation  (9)  that  6  is  determined  by  the  zeros  of  the  cubic  monic  polynomial 


m  =  t^- 


fe  +  6  ,2  .8  +  45^^ 

4  16 


Re’a 

~l6~  ■ 


Assuming  that  Re  a  and  Im  a  are  different  from  zero  we  readily  see  that 


m  =  - 


Re’a 

16 


<0 


/(!)  = 


Im’o 

16 


>0 


so  that  /(f)  has  at  least  one  zero  in  the  interval  0  <  f  <  1.  In  fact,  /(f)  has  three  distinct 
zeros  in  the  interval  0  <  f  <  1.  To  see  this,  let  si,S2,43,«4  be  the  distinct  roots  of  equation 
(5)  which  we  may  write 


[(s  -  Si)(s  -  S2)][(s  -  S3)(3  -  34 )]  =  0 

[(3  -  3i)(3  -  33)][(3  -  32)(3  -  34)]  =  0 

[(3  -  3i)(3  -  34)][(3  -  32)(3  -  33)]  =  0 


Hence  we  may  factor  equation  (5)  into  two  quadratics  in  three  different  ways,  so  that  in 
equation  (6)  we  would  expect  three  distinct  triples  (uj,U2,0).  Further,  since  the  cubic 
equation  /(f)  =  0  has  three  distinct  real  roots,  its  discriminant  is  positive  and  we  can  easily 
write  formulas  for  these  three  roots  (see  for  example  [4,  p.  125]).  Thus  if  we  let  T  be  any 
zero  of  /(f)  we  have  from  equations  (10)  and  (11) 


1  /Re  a  Im  a  \ 

“  2  V  W  VT^) 

_  1  fRe  a  Ima^ 

'‘’"JVvf  ”  VI -t) 


(12) 

(13) 
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where  0  <  T  <  1.  In  particular,  setting 


ca  =  - 


6  +  6  8  +  46+|aP 

-  = - 16 - 


Co  =  - 


Re^a 


4  '  •  16  ’  16 

p  =  cx  -  <^/3  ,  q  =  {24-  9c, Ca  +  27co)/27 


the  zeros  of  /(t)  are 


1  -1  /3g 

u;  =  -  cos  j  ~  4/  — 

3  \2p  V  P  > 


and  we  may  choose  T  to  be  any  one  of  the  latter  expressions. 
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3.  The  Roots  of  the  Quartic 

Noting  equation  (7)  we  write  the  four  roots  of  equation  (6): 

51  =  ^  ^—Hi  cos  6  —  ^4  —  /if  sin  ^  sin  6  +  \Ja  —  nl  cos  6^ 

52  =  ^  cos  6  +  y/4  —  fix  sin  ^  sin  6  —  yj4  —  cos  d'j 

53  =  ^  ^—fi2  cos  6  +  \Ja  —  fil  sin  ^  sin  6  +  \/a  —  fil  cos  0^ 

54  =  ^  ^~/*2  cos  0  —  ^Ja  —  sin  ^  ^Ms  sin  0  —  \Ja  —  /i|  cos  ; 

and  now  applying  the  rotation  equation  (4)  we  obtain  the  four  distinct  roots  of  equation  (1): 


2*  =  exp  i(arg  Sfc  -  ^  arg  a)j  ,  (fc  =  l,2,3,4) 


where 


and 


arg  a  =  tan' 


arg  Si  =  tan 


.j  /  sin  2<^  —  ^1^2  sin  <i>  \ 
\cos  24>  —  ki  ^2  cos  <f> ) 

—fti  v^l  —  T  +  yjA  —  ii\  y/r 


arg  S2  =  tan 


—  fan  ^ 


—  fan  ^ 


^3  =  tan 


arg  S4  =  taa  ' 


-Ml 

Vf 

— 

A 

-M? 

y/i 

-T 

-Ml 

VT- 

-T 

'  — 

A^ 

-Ml 

y/f 

-Ml 

y/f 

+ 

A 

-M? 

yr 

-T 

M2 

f 

+  1 

A^ 

yff 

-M2 

y/f 

+ 

A 

-mI 

VT 

-T 

M2 

VT^ 

T 

—  - 

yPf 

-M2 

y/f 

— 

A 

-Mi 

vT 

-T 

(14) 

(15) 

(16.1) 

(16.2) 

(16.3) 

(16.4) 


In  equations  (16)  fii  and  ft2  are  given  respectively  by  equations  (12)  and  (13)  and  T 
(0  <  T  <  1)  is  any  zero  of  the  cubic  polynomial  /(t).  We  remark  that  the  signs  of  the 
numerator  and  denominator  of  the  argument  of  the  inverse  tangent  function  in  equations 
(15)  and  (16)  must  not  be  disturbed  by,  for  example,  simplification  of  these  equations. 
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4.  The  Grazing  Angle 

From  equation  (14)  we  now  find  the  four  angles: 


V»i  =  2  (“8 
V'2  =  ^  (arg 

V'3  =  ^  (arg  S3 
=  \  (“g 


1  V 

2  "g  “) 

1  . 
2  "g  “) 

1  V 

2  “g  “) 

1  X 

2  "g  «) 


and  the  unique  value  of  that  satisfies  equation  (2)  is  the  grazing  angle  ip.  We  observe 
that  since  ip  must  be  in  the  interval  (0, 7r/2),  we  can  reject  immediately  any  of  the  ipi  not  in 
this  interval. 

It  is  interesting  to  note  that  if  we  choose  the  zero  of  f[t)  to  be 


T  =  2 


extensive  numerical  computations  indicate  that  the  grazing  angle  is  given  by 


^  _  I  5(arg  Si  -  I  arg  q),  Q  <  <p  <  7r/2 
{  |(arg  sj  -  1  arg  a),  7r/2  <  <p  <  rr  , 

Although  the  authors  cannot  prove  this  conjecture  analytically,  it  is  certainly  true;  and 
applied  workers  should  not  hesitate  to  use  it,  since  it  is  a  simple  matter  to  check  numerically 
by  substitution  whether  ip  is  in  fact  the  unique  root  of  equation  (2). 
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5.  Special  Cases 
Since 


we  have 


Re  a  =  -  co6(arg  (3-  ^  arg  a) 

Im  a  =  —  sin(arg  —  -  arg  a) 
a  2 


3  * 

Re  a  Im  a  =  —  sm(2  arg  —  arg  a) 

Q 


so  that  Re  a  or  Im  a  vanish  provided  that 


2  arg  /3  =  arg  a  +  nw  ,  (n  =  0,  ±1,  ±2, . . .) 


Taking  the  tzmgent  of  both  sides  of  this  equation  gives 

Ima_  Re/3Im/3 
Re  a  Re^  ft  —  Im*  /3 

which  yields  after  a  little  computation  the  result  that  Re  a  or  Im  a  vanish  provided  that 
ki ,  fcj,  </)  satisfy  the  equation: 


For  0  <  fci ,  ^2  <  1)  ^  ^2  the  right  side  of  equation  (17)  is  boimded  by  one  so  that  4>  always 

exists  in  the  interval  (0, 7r/2). 

Numerical  computations  show  that  when  equation  (17)  holds,  then  Re  a  =  0  and 
Im  a  <  0;  it  appears  then  that  Im  a  never  vanishes.  Nevertheless,  for  completeness,  we  shall 
discuss  below  the  two  special  cases  Re  a  =  0  and  Im  o  =  0,  since  we  have  not  proved  the 
latter  assertion  but  only  indicated  that  its  truth  is  supported  by  numerical  evidence. 

In  the  case  when  a  is  real,  then  it  is  easy  to  show  from  equations  (6)- (9)  that  the  roots 
of  equation  (5)  are  given  by 

^(-p,  +t  v^4-/i?)  , 
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where 


Ml  =  +  s/a?  -4&  +  8) 

M2  =  -^{a  -  -  46  +  8)  ; 

and  m  the  case  a  (=  ic,  c  real)  is  purely  imaginary  the  roots  of  equation  (5)  are 

^(-Mit+  ^4  -^?)  , 

-  M?)  , 

^(-M2t+  \/4  -  p|)  , 

^(-M2t-  \/4  -  /i|) 


(18.1) 

(18.2) 


where  now 


Ml  =  ^(c+  y/<?  +  46  +  8) 


M2  =  2^^”  Vc*  +  46  +  8)  . 


(19.1) 

(19.2) 


Therefore,  if  0  <  fcj,  ^2  <  1,  fci  ^  fc2,  0  <  0  <  ir  are  such  that  equation  (17)  holds,  then 
in  the  case  a  is  real  we  have 


^1  =  ^ 


tan-  (  V'  “  "‘'l 

1 

-  O  arg  O 

\  -Ml  ; 

2 

tan  ^  [ 

v4-Mr 

\  1 

\ 

/  2  J 

tan“* 

/  v/^-Ml'j 

1 

-  -  arg  Q 

- 

V  / 

2 

tan“^  j 

-  v/4^' 

\  1 

-  arg  Q 

\ 

-M2  ; 

1  2 
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where  arga  and  are  given  respectively  by  equations  (15)  and  (18);  and  in  the  case 

where  a  =  ic  (c  real)  we  have 


V'l  =  - 

2 


tan 


-1  /  -Pi 


V’2  =  - 

2 


va  2 


tcLn 


-1  /  -Pi 


P?> 


tan 


-v/r: 

-1  /  -P2 


P?. 


tan 


-1  /  -P2 


P2> 


v/4^ 

where  fii  and  /xj  axe  given  by  equations  (19). 


P2> 


-  larga 


^arga 


^arga 


-larga 
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6.  Conclusions 

In  the  present  paper  we  have  shown  that  it  is  possible  to  exhibit  the  grazing  angle  of 
specular  reflection  as  one  of  four  angles  derived  from  the  zeros  of  a  certain  self-inversive 
quartic  polynomial  whose  coefflcients  are  complex.  Thus  the  grazing  angle  may  now  be 
found  deterministically,  whereas  heretofore  it  had  to  be  computed  either  by  approximations 
or  numerical  procedures. 

The  preceeding  analysis  shows  further  that  factoring,  in  a  useful  way,  a  quartic  polynomial 
whose  complex  coefflcients  depend  on  even  a  small  number  of  parameters,  is  not  always  a 
trivial  matter,  even  though  the  methods  of  Cardan  and  Ferrari  have  been  known  for  centuries. 
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